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Here we evaluate the many-body entanglement properties of a generalized SU(n) valence bond 
solid state on a chain. Our results follow from a derivation of the transfer matrix of the system which, 
in combination with symmetry properties, allows for a new, elegant and straightforward evaluation 
of different entanglement measures. In particular, the geometric entanglement per block, correlation 
length, von Neumann and Renyi entropies of a block, localizable entanglement and entanglement 
length are obtained in a very simple way. All our results are in agreement with previous derivations 
for the SU(2) case. 
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Introduction.- The study of entanglement in strongly 
correlated systems has proven fruitful to understand new 
phases of quantum matter and new types of quantum or- 
der [1]. In this respect, there has been growing interest 
in quantifying entanglement in the ground state of quan- 
tum many-body systems in one spatial dimension [2] . At 
criticality, the entanglement in these systems diverges, 
in turn obeying precise scaling laws orchestrated by the 
underlying conformal symmetry. Away from criticality, 
though, the existence of a finite correlation length and 
a non-zero gap to excitations forces entanglement to re- 
main finite. 

The archetypical example of a quantum spin chain with 
a gap is the spin-1 AKLT model, introduced in Ref. [3] by 
Affleck, Kennedy, Lieb and Tasaki. This model is invari- 
ant under rotations, that is, SU(2) operations. Moreover, 
its ground state is a valence bond solid (VBS) that ad- 
mits a representation in terms of a Matrix Product State 
(MPS) [ ], and is closely related to the Laughlin state [5] 
and the fractional quantum Hall effect [ ]. This scenario 
has been recently generalized to other symmetry groups 
such as SO(n), SU(n) and Sp(2n) [7-12]. As for the be- 
haviour of entanglement in these generalizatons, not too 
much is known. Derivations have been carried out for 
the correlation length [8, 9] as well as von Neumann and 
Renyi entropies [10, 11] of SU(n) VBS states on a chain, 
but these involve a number of technicalities that make 
them quite lengthy. 

In this paper we provide an elegant and straightfor- 
ward evaluation of the many-body entanglement prop- 
erties of the above SU(n) valence bond solid state on 
a chain. In particular, we derive unknown quantities 
such as the geometric entanglement per block [13, 14], 
but also re-derive other quantities such as the correla- 
tion length, von Neumann and Renyi entropies of a block 
in a significantly simpler way than previous derivations 
[10, 11]. Our calculations are novel in many aspects for 
SU(n) VBS states and are based on a proper understand- 
ing of (i) the structure of the MPS transfer matrix of a 
block, and (ii) the constraints imposed by SU(n) sym- 
metry. Furthermore, we also consider the localizable en- 
tanglement [15] and prove that the entanglement length 



diverges in the thermodynamic limit [16]. 

The SU(n) VBS state: We start by reviewing the con- 
struction of the one-dimensional VBS state with SU(n) 
symmetry from Refs.[10, 11]. This is the SU(n) analogue 
of the ground state of the spin-1 AKLT chain, and can 
be represented for N sites as 

I*) = (S^WS 1 ) \^of' 1 ¥oo) [l - 2] ---\^o) [N - N+1] ■ 

In the above equation, |<I ) oo)' r:r+1 ' is a singlet state of two 
'virtual' particles at sites f and r + 1, one in the conju- 
gate (□) and one in the fundamental (□) representations 
of SU(n), both of dimension n. Operator wj[^ is the 
projector onto the (n 2 — l)-dimensional adjoint subspace 
of □ <8> □. This adjoint subspace is the physical space at 
site r. As explained in Refs.[9-ll], state |^) is the ground 
state of a gapped parent Hamiltonian with SU(n) symme- 
try. In fact, for n > 2 state |\&) violates parity symmetry 
since the fundamental and conjugate representations of 
SU(n) are different. This implies that for open bound- 
ary conditions it is the unique ground state of its parent 
Hamiltonian, whereas for periodic boundary conditions 
there are two degenerate ground states (one for each par- 
ity sector) depending on whether we choose the singlets 
|$ 00 ) for the tensor products □ (g) □ or □ ® □ . The VBS 
state I*) is represented diagramatically in Fig.(l.a). 

Let us elucidate further the structure of the above 
state. First, notice that by assigning a basis \j) to 
the fundamental representation and \j) to the conju- 
gate (with j = 0,1, ...,n— 1), the singlet |*l>oo) can 
always be represented as a maximally entangled state, 
i.e. |$oo) = rfe Sj=o li)U)- Second, we remind that the 
tensor product of the fundamental and conjugate rep- 
resentations of SU(n) obeys □ ® □ = singlet adjoint, 
where by 'singlet' we mean the singlet representation (of 
dimension 1) and by 'adjoint' we mean the adjoint repre- 
sentation (of dimension n 2 — 1). In terms of a complete- 
ness relation for vector spaces, this implies that 

I„ «) I„ = |$00><$00|+W„* , (2) 

where I„ = J3?=o b')01 i s ^ ne 71 x 71 identity matrix be- 
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FIG. 1: (a) Structure of the SU(n) VBS state |*) for N sites 
and open boundary conditions, (b) MPS structure of the state 
(up to an overall normalization constant). At the boundaries 
there are SU(n) degrees of freedom a, b = 0, 1, . . . , n — 1. (c) 
MPS transfer matrix A(l). (d) Transfer matrix of a block 
A(L) = A(1) L . 



tween two fundamental spaces, and analogously for E„. 
Thus, we can understand operator W a gj — I„ ® I n — 
|$oo)(^ool as a projector that removes the singlet com- 
ponent of the tensor product □ <g> □. 

At this point, we wish to introduce a very convenient 
notation in terms of tensor network diagrams: we will 
draw a line every time we wish to represent an identity 
matrix. This identity matrix will be between two funda- 
mental spaces, two conjugate spaces, or one fundamental 
space and one conjugate space, with no distinction from 
the diagrammatic point of view. Using this notation, 
Eq.(2) reads 
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W ad 



(3) 



In the above expression, the left hand side is nothing but 
a diagram representing I„(g)I„. To clarify the meaning of 
the right hand side, notice that the singlet state must be 
represented as |$oo) = [l/v 7 ^ We thus represent 

the projector in the singlet subspace as — , which is 
the first term in the right hand side of Eq.(3). With this 
notation it is also possible to represent the MPS struc- 
ture of state | in a clear and elegant way as shown in 
Fig.(l.b), where the MPS matrices are just proportional 
to the projectors W a dj at every site [22]. 

Formally, we can think of Eq.(3) as the completeness 
relation for the tensor product of 'left' and 'right' SU(n) 
representations. From a geometrical perspective, though, 
one can also have the tensor product of 'up' and 'down' 
representations. In this case, one has the relation 
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(4) 



with W^ d j the projector onto the adjoint subspace of the 
considered tensor product. Both Eqs.(3) and (4) will turn 
out to play a key role in our derivations. 

Transfer matrix of a block: All the calculations in 
this paper rely on the derivation of the MPS transfer 



matrix A(i) for state \^f). This is defined in terms 
of the MPS matrices W adj as A(l) = {l/n)W ad3 W' c 



see Fig.(l.c). Since W ad j is a projector, we have that 
w ad 3 Wl dj = W adj W ad:j = W 2 dj = W ad] , which can be 
checked easily using e.g. Eq.(3). So, in the end, we ob- 
tain the remarkably simple expression 
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We now wish to compute the transfer matrix A(L) of 
a block of length L. As shown in Fig.(l.d) this is given 
by A(L) — A(1) L , and thus can be obtained easily af- 
ter diagonalizing A(l) with respect to 'left' and 'right' 
spaces. This diagonalization is indeed straightforward if 
one realizes that the transfer matrix A(l) in Eq.(5) can 
be rewritten using Eq.(4) as 



A(l) = 



1 - 



1 



-1 



v f adj 



(6) 



which is a sum of orthogonal projectors in singlet and 
adjoint subspaces, and is thus diagonal. Its largest eigen- 
value is non-degerenate and is given by Ai = (1 — 1/n 2 ) 
corresponding to the singlet eigenspace. There is also a 
second eigenvalue A2 = — 1/n 2 with degeneracy n 2 — 1 
corresponding to the adjoint eigenspace. Using this, the 
transfer matrix of a block reads 



A(L) 



1 



1 L 



-1 



' y adj 



(7) 



As expected, the above expression for A(L) is man- 
ifestly SU(n) invariant. To see this, just notice that 
thanks to Eq.(4) one can regard A(L) as a combination 
of tensor products ofnxn identity operators, which is of 
course invariant under any SU(ro) transformation. This, 
in fact, is a consequence of the SU(n) invariance of the 
quantum state |\E f ). 

Correlation length and state norm. - Some properties of 
the SU(n) VBS state |*) can already be computed just by 
looking at Eqs.(6) and (7). For instance, the correlation 
length is given by £ c = -1/log |A 2 /Ai| = l/log(n 2 - 1), 
in agreement with results present in the literature [8, 9]. 
Furthermore, we can also compute the norm ( V E'| , I / ) of the 
state, which will be needed in the forthcoming calcula- 
tions. For a system of size N and open boundary condi- 
tions this is given by (^1^) = nX^ , whereas for periodic 
boundary conditions (*|*) = Af + (n 2 - 1)A^, with A1 
and A2 defined as above. In the thermodynamic limit the 
norm scales as A^ -1 log (^l^) ~ log Ai = log (l — 1/n 2 ) 
regardless of the boundary condition. 

Geometric entanglement per block.- We now con- 
sider the behaviour of the geometric entanglement per 
block of length L, which we call £{L) [ ]. This 
is defined in the thermodynamic limit as £{L) = 
lim (JV/i) ^ 00 (7V/L)- 1 ^(*), with E(9) = -log(|A max | 2 ) 



where the maxi- 



and |A max | 2 = max (<S>\V) / y/(*\*j 

mization is done over all possible product states |$) of 
contiguous blocks. This measure of entanglement quanti- 
fies the resemblance of the quantum state 14") to a prod- 
uct state of the blocks. Thus, studying how £ (L) scales 
with L allows us to understand how similar \^>) is to 
a product state under successive coarse-grainings of the 
system [14]. In our case, since the system is gapped, we 
expect £ (L) to saturate in some constant for L ^> £c- 
This will turn out to be the case, as we shall see shortly. 

Let us now perform the calculation of £{L). Our 
derivation follows similar lines as the calculations in 
Ref.[17]. First of all, since the system is antiferromag- 
netic the basic unit cell is 2, and thus we consider here 
the case of blocks of even length L only [23]. Sec- 
ond, as explained in Refs.[17] and [18], for an MPS 
with diagonalizable transfer matrix A(l) the quantity 
|A max | 2 is given in the limit TV > 1 by |A max | 2 ~ 
(max\(r\*(r\A(L)\r)\r)*\) N/L where |r), \r)* , (r\* 
and (r| are normalized vectors in the four spaces con- 
nected by A(L) (respectively right-up, right-down, left- 
up and left-down), and the maximization is done over 
these vectors. Generally speaking, carrying this maxi- 
mization is cumbersome. However, an important sim- 
plification occurs in our case since A{L) is just a com- 
bination of tensor products of identity operators. This 
implies that the maximization does not depend on vector 
\r) at all, which is a consequence of the SU(n) symmetry 
of A(L) [24]. Using Eq.(7) and our previous result for 
the norm, one can then see that 



, ,'(1 - l/n 2 ) L (l/n) + (-l/n 2 ) L (l - l/n) 
(1- l/n 2 ) L 



N/L 



(8) 

which means that the geometric entanglement per block 
is ultimately given by 

£{L) = logn- log (l + (n- l)e~ L/Cc ) L even , (9) 

where we have used our previous result for the correlation 
length £(7. As expected, the above formula matches the 
SU(2) case calculated in Ref.[17] for n = 2. Also, for 
L ^> £c the geometric entanglement saturates in £(L ^> 
1) ~ logn, as expected for a gapped system. We also 
see that for n ^> 1 the geometric entanglement increases 
linearly with the size of the block, £ (L) ~ L/£,c, before 
reaching the saturation regime for large L. 

Von Neumann and Renyi entropies of a block.- Our 
derivation of the transfer matrix of a block allows us 
to compute as well the Von Neumann and Renyi en- 
tropies of a block, S(L) and S a (L). These are defined 
respectively as S(L) = — J2 a ^'a^°E^'a an d S a (L) 
: 1 — a) -1 log (X] a (Aa) a ) for a 7^ 1 and a > 0, and where 
X' a are the eigenvalues of the reduced density matrix of a 
block of size L, which we call pl- The set of Renyi en- 




FIG. 2: Diagrammatic representation of several reduced den- 
sity matrices sharing the same spectrum (up to normalization 
constants), (a) Of a block of length L, called pL- (b) Of the 
environment of a block of length L for the environment as in 
(a), (c) Of a block of length L once the environment is com- 
puted (dominant left and right eigenvectors of A(l), which 
are singlets), (d) Of the environment of a block of length L 
for the environment as in (c), equal to A(L). 



tropies contains lots of information about the entangle- 
ment in the system. Both the von Neumann and Renyi 
entropies were already computed in Ref.flO] for the SU(n) 
VBS state I*) using a different method. Here, though, 
we re-derive them by using the formalism introduced ear- 
lier, which allows for a new, simple and elegant derivation 
in just a few lines. 

The calculation of these entropies amounts to being 
able to diagonalize the reduced density matrix of a block 
Pl- In terms of tensor network diagrams, this reduced 
density matrix is represented in Fig. (2. a, c) (up to nor- 
malization). The key observation is that the eigenval- 
ues oi pl must be equal to those of the reduced density 
matrix of the environment of the block, see Fig.(2.b,d). 
Thus, up to a normalization constant, we can write 



PL A(L) 



(10) 



where «-» means that they share the same eigenvalues if 
we diagonalize A{L) with respect to the 'up' and 'down' 
spaces. This diagonalization can be done by using again 
the decomposition given in Eq.(3). We obtain the diago- 
nal expression 
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1 - 
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1 

n 
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W ad3 . (11) 



Defining p n (L) = (-l/(n 2 - 1)) L = (-1) L exp (-L/fc), 
we find the normalized eigenvalues 



K = - 2 {l + (n 2 -l)p n (L)) 



An = 



r(l-Pn(i)) , 



(12) 



where X[ has degeneracy 1 (singlet eigenspace) and X' 2 
has degeneracy n — 1 (adjoint eigenspace). These eigen- 
values coincide exactly with the ones given in Eq.(24) 
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of Ref.[10], and from them the von Neumann and Renyi 
entropies follow inmediately. Notice that this derivation 
is also far more straightforward than the one in Ref.[10] 
since, essentially, we just used repeated applications of 
Eq.(2) in different contexts. As expected, for the n = 2 
case one recovers the expressions for the spin-1 AKLT 
state in Ref.[19]. Furthermore, in the limit L ^> £c we 
see that ~ \' 2 ~ 1/n 2 . In this limit the reduced den- 
sity matrix p L is thus proportional to some projector on 
an Tridimensional subspace, and the von Neumann and 
Renyi entropies saturate as S(L) ~ S a (L) ~ 21ogn, as 
was also observed in Ref.[10]. Interestingly, this satu- 
ration value is twice the one that we obtained for the 
geometric entanglement per block in the same limit. 

Localizable entanglement and entanglement length.- As 
in the SU(2) case [ ], it is also possible to compute 
exactly the localizable entanglement and entanglement 
length of this system. By definition, the localizable en- 
tanglement is the maximum entanglement that can be lo- 
calized, on average, between two parties of a multipartite 
system, by doing local measurements on the other par- 
ties. This brings to the concept of entanglement length 
£e as the length scale at which it is possible to create a 
maximally entangled state between two parties by doing 
measurements on the rest of the parties. In our case the 
parties are the different sites of the system, correspond- 
ing to adjoint representations of SU(n). 

Our derivation of £e is a generalization of the one in 
Ref.[16]. Specifically, we consider the generalized Bell 
basis for two qu-nits |$, p ) = -J= Y^ 1 e 27r *>|j © p )|J) j 
with l,p — 0,1, ... ,n — 1. Notice that for I = 0,p = we 
recover the singlet state |$oo)- A measurement in this 
basis achieves perfect teleportation for the state of one 
qu-nit [20]. The key observation is that, of the above 
states, |$oo) is the singlet, whereas the rest of states are 
an orthonormal basis of the adjoint (n 2 — l)-dimensional 
subspace. Thus, a measurement at each site on this basis 
for l,p ^ amounts to a generalized Bell measurement 
on the virtual particles. As in Ref.[16] these measure- 
ments implement entanglement swapping, which means 
that the localizable entanglement between any pair of 
sites is always maximum. In particular, such measure- 
ments allow to create a maximally entangled state be- 
tween the two qu-nits at the boundaries of the system. 
Therefore, for N sites we have that £e ~ N, which di- 
verges in the thermodynamic limit. Importantly, and as 
in Ref.[I6], the entanglement length £g diverges while 
the correlation length £c is finite. 

Conclusions and outlook.- Here we have shown how to 
obtain the entanglement of a VBS state on a chain with 
SU(n) symmetry in a novel, elegant and straightforward 
way, much simpler than previous derivations. We believe 
that similar techniques to the ones in this paper could 
also be useful in the study of one-dimensional VBS states 
with other symmetries such as SO(n) and Sp(2n) [7, 12], 
as well as SU(n) VBS states in higher dimensions [ ]. 
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